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We present a measurement of the branching fraction for the decay of the neutral B meson into 
the final state J/if)-K + ir~ . The data set contains approximately 56 million BB pairs produced at 
the T(4S) resonance and recorded with the BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage ring. The result of this analysis is B(B° — > J/ip-K + ir~) = (4.6 ± 0.7 ± 0.6)xl0~ 5 , where 
the first error is statistical and the second is systematic. In addition we measure B(B° — > J/ip p°) = 
(1.6±0.6 ±0.4) x 10~ 5 . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



In the Standard Model, the decay B a —> J/tp p° 
can give rise to CP-violating asymmetries (directly and 
through B°-B° mixing). Therefore it is interesting to 
study the decay mode _B° — > J/ipTT + n~ to understand the 



J/ipp component in the final state. Since these decays 
are Cabibbo and color suppressed, they are a sensitive 
probe of rare and exotic physics processes, such as pen- 
guin contributions and box diagrams containing charged 
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Higgs bosons. These effects may appear as deviations of 
the branching fraction from the Standard Model predic- 
tion of B(B° -> J/tPtt+tt-) = (4.6 ±0.8) x 1(T 5 §. This 
decay mode has not previously been observed. CLEO 
quotes an upper limit of B(B° — > J/ipp°) < 2.5 x 10~ 4 
at the 90% confidence level Q . Here we present the first 
measurement of B(B° — » J/07r + 7r _ ). 

The data used in the present analysis were collected 
at the PEP-II storage ring with the BABAR detector, 
described in detail elsewhere ||. Charged particles 
are detected, and their momenta measured, with a 40- 
layer drift chamber (DCH) and a five-layer silicon vertex 
tracker (SVT), both operating in a 1.5 T solenoidal mag- 
netic field. Surrounding the DCH is a detector of inter- 
nally reflected Cherenkov radiation (DIRC), and outside 
this is a CsI(Tl) electromagnetic calorimeter (EMC). The 
iron flux return of the solenoid is instrumented with re- 
sistive plate chambers (IFR). The data sample used for 
the analysis contains approximately 56 million BB pairs, 
corresponding to a luminosity of 51.7fb _1 recorded near 
the T(4S) resonance. An additional 6.4fb _1 , recorded 
approximately 40 MeV below the T(4S) peak, were used 
to study continuum backgrounds. 

Events containing BB pairs are selected based on track 
multiplicity and event topology [Q . At least three tracks 
are required to originate near the nominal beam spot, 
with polar angle in the range 0.41 < #i a b < 2.54 rad, 
transverse momentum greater than 100 MeV/c, and a 
minimum number of DCH hits used in the track fit. To 
reduce continuum background the ratio of second to ze- 
roth Fox- Wolfram moment, i?2 = H 2 /H , is required to 
be less than 0.5. The sum of charged and neutral energy 
must be greater than 4.5 GeV in the laboratory frame. 
The primary vertex of the event must be within 0.5 cm 
of the average measured position of the interaction point 
in the plane transverse to the beamline. 

The J/tjj is reconstructed in the e + e~ and final 
states. Electron candidates must satisfy the requirement 
that the ratio of calorimeter energy to track momentum 
lies in the range 0.75 < E/p < 1.3, the cluster shape and 
size are consistent with an electromagnetic shower, and 
the energy loss in the DCH is consistent with that for an 
electron. If an EMC cluster close to the electron track is 
consistent with originating from a bremsstrahlung pho- 
ton, it is combined with the electron candidate. 

Muon candidates must satisfy requirements on the 
number of interaction lengths of IFR iron penetrated 
(N\ > 2), the difference between the measured and ex- 
pected interaction lengths penetrated (\N\ — N^ xp \ < 2), 
the position match between the extrapolated DCH track 
and the IFR hits, and the average and spread of the num- 
ber of IFR strips hit per layer. 

Pion candidates are accepted if they originate from 
close to the beam spot and are not consistent with being 
a kaon. The algorithm uses dE/dx information from the 
SVT and DCH, and the Cherenkov angle and number of 



photons from the DIRC. 

Tracks are required to lie in polar-angle ranges where 
particle identification efficiency is measured with known 
control samples. The allowed ranges are 0.41 < #i a b < 
2.41 rad for electrons, 0.30 < 6i a b < 2.70 rad for muons, 
and 0.35 < 9\ a b < 2.50 rad for pions, which corre- 
spond approximately to the geometrical acceptances of 
the EMC, IFR, and DIRC, respectively. 

Identified electron and muon pairs are fit to a common 
vertex and must lie in the J/tp invariant mass interval 
2.95 (3.06) to 3.14GeV/c 2 for the e+e' (/i + M~) channel. 

B candidates are formed by combining a J/ip candi- 
date with a pair of oppositely-charged pion candidates 
consistent with coming from a common decay point. We 
also require the positions of the vertices of the lepton 
pair and the pion pair to be consistent. Further selec- 
tion requirements are made using two kinematic vari- 
ables: the difference, AE, between the energy of the can- 
didate and the beam energy £^" am in the center-of-mass 
frame, and the beam-energy substituted mass, toes = 
\A^bcam) 2 — (p'b 1 ) 2 - After applying the loose require- 
ments 5.2 < m ES < 5.3GeV/c 2 and \AE\ < 0.12 GeV, 
approximately one-quarter of events contain more than 
one B° candidate, from which we keep the one with the 
smallest \AE\. The distribution of the candidates in AE 
and toes is shown in Fig. |l|. For the final signal sam- 
ple, we require |m ES - 5279.0 MeV/c 2 < 9.9 MeV/ c 2 and 
\AE\ < 39 MeV, which correspond to 4er and 3tr ranges 
in the resolutions for toes and AE. After all selection 
criteria have been applied, 213 events remain. 



0.1 
0.05 


-0.05 
-0.1 



<• .-A . : -ft :>w> -^M-r ■ v; 




FIG. 1: Signal for B° — » J/ip-K + -K~ . The upper plot shows the 
distribution of events in the A_E-mES plane, where the box 
represents the final selection criteria. The lower plot shows 
the distribution in mES of events with \AE\ < 39 MeV, where 
the dashed (solid) line corresponds to events in the K% (non- 
K° s ) region in M(-k + -k~) (0.45-0.55 GeV/c 2 ). The vertical lines 
represent the final selection. 



An unbinned, extended maximum-likelihood Q fit is 
performed on the invariant mass distribution of the two 
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pions for the selected events, to determine the vari- 
ous contributions to the B° — > J/07r + 7r _ events. We 
consider five categories: (i) B° — » J/tpp° events; (ii) 
B° -> J/iPK°(K° -> TT+TT") events; (hi) 5° -> J/^tt+tt- 
(non-/o° signal) events; (iv) background from events with- 
out a real J/0 ; (v) inclusive- J/0 background from events 
containing a real J/0 . A probability density function 
(PDF) is constructed for each of these five cases. The 
total PDF is then formed from the sum of the five PDFs 
and fit to the data. The B° — > J/tp Kg mode is not con- 
sidered to be signal for the purposes of determining the 
branching fraction for B° — > J/07r + 7r~ . 

The PDF used to model the 73° — > J/ipp° mode is a 
relativistic P-wave Breit-Wigner function 0: 
F p (m) = (mr(m)P 2i * ff+1 )/((m^-m 2 ) 2 + TOp"(m) 2 ), 

where F(m) = r (^) 3 (=f)(i±ff|). g(m) is the pion 
momentum in the di-pion rest frame, with qo = q(m p ). 
m = M(7T + 7r - ) is the two-pion invariant mass and P 
is the J/0 momentum in the -B rest frame. m p — 
770MeV/c 2 , T = 150MeV/c 2 , and m n = 140MeV/c 2 . 
P e ff is the effective orbital angular momentum between 
the J/0 and the p°, which can take any value between 
and 2 and so is allowed to float in the fit. R is the 
Blatt-Weisskopf barrier- factor radius . The fit is per- 
formed with R equal to two values (0.5 and 1.0 fm) and 
the results of the two fits are averaged. 

The PDF for the B° -> J/tpK° s mode is a single Gaus- 
sian function with the mass and width fixed to values 
obtained by fitting a sample of simulated J/tpK® events. 
Allowing these parameters to vary in the final M(it + it~) 
fit does not change the results. 

The PDF used to model the B° -> J/tpir+Tt~ (non- 
p° signal) contains a three-body phase space factor 
q(m)P(m) and a factor of P(m) 2 motivated by angular 
momentum conservation: Fp^(m) — q(m)P(m) 3 . If the 
7r + 7r~ is in an S'-wave, angular momentum conservation 
results in a factor of P(m) 2 , while a D-wave yields the 
second power of P(m) or higher. 

The PDF for the Af(7r+7r~) distribution for back- 
ground events without a real J/0 is derived from a fake- 
J/0 sample selected in data as described above except 
that at least one of the lepton candidates must fail the ap- 
propriate particle identification requirements. A Monte 
Carlo study confirms that the M(tt + tt~) distribution ob- 
tained with this procedure correctly describes the shape 
of the non-J/0 background. The resulting distribution is 
parametrized using the sum of two Weibull functions 
and a Breit-Wigner. The Breit-Wigner describes the p° 
component of the non- J/0 background. 

The PDF for the M(7r + 7r~) shape for background 
events containing a real J/0 is obtained from a simulated 
B — ► J/tpX sample equivalent to a luminosity of 81 fb -1 . 
Events in which the system X is 7r + 7r~ (non-resonant), 
p°, or ifj (ir + ir~) are removed from the sample. The 
resulting shape is described by a Weibull function. 



The normalization of the background components is 
obtained from samples in data and simulation. The level 
of non- J/0 background is obtained from sidebands of the 
J/0 mass distribution in data. The rriES distribution 
for these sideband candidates is then fit to an ARGUS 
function [|| to determine how many events pass the final 
selection criterion. Scaling to the equivalent background 
in the J/0 mass region, using an exponential to describe 
the background shape in the J/0 mass distribution, the 
expected non-J/0 background is found to be 35.7±1.2 
events. 

The level of inclusive- J/tp background is obtained from 
the distribution of to_es f° r events in the AE signal region 
in both data and simulation. In each case the rriES distri- 
bution is parametrized by a Gaussian function (to repre- 
sent signal or peaking background) and an ARGUS func- 
tion. Peaking background originates from B — * J/tpX 
decays such as B — > J/tpK*, B + — > J/0 p + , and B — » 
J/i/jKi, that accumulate near tties = 5.279 GeV/c 2 . 

The non-peaking component of the inclusive- J/tp back- 
ground is determined by subtracting the non- J/tp contri- 
bution, on the basis of the scaled sideband events de- 
scribed above, from the total ARGUS background in 
data. The peaking component is determined from the 
Gaussian part of the m es distribution in B — -» J/tp X sim- 
ulation, where events with X = tt + it~ (non-resonant), p° 
and Ks(tt + tt~) have been removed. The sum of peaking 
and non-peaking components of the inclusive- J/tp back- 
ground is found to be 61±11 events, of which the peaking 
component comprises 6 events. Thus any associated un- 
certainties, such as branching fractions used in the J/tpX 
simulation, will not contribute significantly to the final 
systematic uncertainty. 

The branching fraction is obtained from 



B{B° -> J/07T+7T- 



N 



J/lp7 



%x fj# „xfi(J/^<+f-) : 



(1) 

where is the total signal yield obtained from 

the fit, N B o is the total number of B° and B° in 
the data sample [El, and £j/^ T7r is the signal efficiency. 
The J/tp branching fraction B(J/tp — > £ + £~) is fixed to 
11.81% pj. We assume that the branching fraction for 
Y(4S) -> B°B° is one-half. 

The signal efficiencies for all requirements apart from 
particle identification criteria are derived from simula- 
tion. Lepton and pion identification efficiencies are de- 
termined with samples of known muons, electrons and 
pions in the data from the following processes: /i + /i~7, 
p+p-e+er, e+e", e+e~7, D*+ -► D°ir+ {D° -> K~tv+), 
and K® —* tt + tt~ . The efficiencies are determined as a 
function of momentum, and polar and azimuthal angle. 
We find e(J/tpp°) = (27.1 ± 0.3)% and e(J/tPir + ir~ , non- 
resonant) = (27.0 ± 0.3)%. The final corrected signal 
efficiency of(27.1±0.2)%is taken as the average of the 
J/tp p and J/07r + 7r~ (non-resonant) efficiencies, where 
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the error is from Monte Carlo statistics. 

A likelihood fit is performed on the M(ir + ir~) distribu- 
tion in data with the normalization of the non- J/ip back- 
ground fixed to 35.7 events and the inclusive- J/ip back- 
ground to 61. Thus only the yields for J/ipp°, J/ipn + ir~ 
(non-p signal), and J/ipK® events are allowed to vary. 
The results of the fit are overlaid on the data points in 
Fig. ||. The goodness-of-fit \ 2 is 33.4 for 38 data points. 





1 






BaBar 











0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

M(7t + 7t-) (GeV/c 2 ) 



FIG. 2: Distribution of the invariant mass M(tt + -k~) for 
events passing all selection criteria. The solid line is the re- 
sult of the unbinned likelihood fit. The dashed line represents 
the sum of background and non-p° signal components. The 
dotted (dot-dashed) line shows the total (inclusive- J/ip ) back- 
ground. The spike corresponds to B° — > J/ip Kg events. 



The result of the fit is 84 ± 13 signal events, of which 
28±10 are in the p° component and 55±15 are in the non- 
p° signal component. The number of events in the 
component is 28 ±5. Inserting the result into Eq. [l] yields 
the branching fraction B(B° -> J/ipn+n-) = (4.6±0.7)x 
10~ 5 , where the error is statistical. 

The signal yield can be checked by counting the 
number of events passing all the selection criteria and 
subtracting the estimated numbers of background and 
J/ipK® events. This method gives 87 ± 15 J/ipir + ir~ 
events, where the error is statistical. 

The systematic errors on the final branching fraction 
measurement arise from uncertainties on the signal effi- 
ciency, fitted yield, number of BB pairs produced, and 
J/ip — » £ + £~ branching fraction. N BB is known to 1.1% 
with the dominant contribution to the uncertainty com- 
ing from the error on the efficiency for the B°B° selec- 
tion. B{J/ip -> £+£~) is known to 1.7% (fractional) @. 

The uncertainty on the pion identification efficiency is 
1.8% per pion. Contributions to this error come from 
the limited size of the data sample used to determine 
the efficiency, the uncertainty on the kaon contamination 
in the sample, and residual differences between the effi- 
ciencies for the known pions in data and for pions from 
B° — > J/i/jtt + it~, determined from Monte Carlo simula- 
tion. 



Uncertainties on electron and muon particle identifi- 
cation efficiencies come from studies using B — > J/ipX 
events in data. Fits to the M(J/ip) distribution in these 
events, under different selection criteria, give estimates 
of the electron and muon identification efficiencies and 
their errors. The overall error for the identification cri- 
teria used in this analysis is 1.3%. 

The uncertainty on the determination of the tracking 
efficiency is 1.3% per track, and is summed for the four 
tracks from the B° decay. The efficiency of the con- 
vergence requirement on the tt + it~ vertex fit has been 
studied with a sample of ip(2S) — ► £ + £~ decays. Data 
and simulation are found to be in good agreement, with 
an associated systematic error of 1%. The unknown p° 
helicity in the J/ip p° component of the final sample in- 
troduces a systematic error on the efficiency of 2.5%, as 
determined from the efficiency variation between simu- 
lated samples with helicity and 1. The limited amount 
of simulated data leads to an uncertainty in signal effi- 
ciency of 0.7%. To determine the effect of the signal and 
background shapes and the background yields on the fit- 
ted yields, the fixed parameters of these PDFs are varied 
within their uncertainties, allowing for correlations. This 
produces a total systematic error due to fit parameter 
variation of 9.7%, which is dominated by the errors in 
the background yields. The final fit neglects resonances 
such as / (980),/ 2 (1270) and p°(1450). Allowing for the 
addition of such terms in the likelihood function results 
in a systematic uncertainty on the yield of 2.1%. The 
total systematic uncertainty from all sources is found to 
be 12.3%. 

The analysis is repeated with variations in the selec- 
tion criteria. Taking into account statistical correlations 
between the results, we find that variations are consis- 
tent with statistical fluctuations due to the addition or 
removal of some of the events in the sample. 

The branching fraction can be measured separately 
for events containing a J/ip — > e + e~ or a J/ip — > 
candidate. The results from these subsamples 
are B(B° J/ijjn+ir-) ee = (5.3 ± 1.1) x 1CT 5 and 
B(B° -> J/ipn+n-)^ = (4.0 ± 1.0) x 10" 5 , where the 
errors are statistical. 

Another way to model the backgrounds is to use a 
smoothing algorithm on the simulated B — > J/ipX and 
fake- J/ip data events, rather than impose a parametriza- 
tion. The resulting PDFs follow fluctuations and check 
how strongly the fitted signal yields depend on the cho- 
sen method of describing the backgrounds. Changing the 
background modeling in this way alters the total fitted 
yield by less than one event. 

The M(ir + Tr~) distribution shows a clear peak at the 
p° mass. The fit result of 28 ± 10 events for the p° sig- 
nal leads to a branching fraction of B(B° — ► J/ipp°) — 
(1.6±0.6 (stat)±0.4 (syst)) x 10~ 5 . The systematic error 
includes a contribution from the effect of using an alter- 
native PDF to describe the non-p° signal. The shape is 
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from a polynomial fit to data recorded in it-it scatter- 
ing experiments jllj] and thus provides an empirically- 
derived shape, in contrast to the default non-p signal 
PDF, which is based on a phase-space assumption. The 
assumption that the non-p° signal is predominantly S- 
wave, and therefore interference with the p° can be ne- 
glected, has been checked on data. A significant S- 
wave contribution means that the leptons from the J/ip 
have a helicity angle distribution oc sin 2 (9j^ ). For 
events in data with M(ir + ir~) > 1.1 GeV/c 2 , we subtract 
the helicity cosine distribution for events with tues < 
5.27MeV/c 2 from the distribution for events in the signal 
hues region and find that the shape of the resulting dis- 
tribution is consistent with sin {6m)- Interference be- 
tween S- and P-wave signal components integrates out 
in the M(tt + ii~) projection, as long as the acceptance 
is symmetric in the cosine of the di-pion helicity angle, 
Q tttt ■ Studies using simulated non-resonant 5- wave events 
show that there is no significant odd component to the 
acceptance function in cos(0 W7r ). Consequently, there is 
no such interference contribution to the ir + ir~ mass dis- 
tribution 

In summary, the branching fraction for B° meson 
decay to the final state <//07r + 7r _ has been measured 
for the first time. The result, B(B° -> J/^Tr+Tr") = 
(4.6 ± 0.7(stat) ± 0.6(syst)) x 10~ 5 , is consistent with 
the Standard Model prediction |Q . In addition, the tech- 
nique of fitting the M(7r + 7r~ ') distribution allows a mea- 
surement of the branching fraction for the J/ip p° compo- 
nent. The result is B(B° -> J/ipp°) = (1.6 ± 0.6 (stat) ± 
0.4(syst)) x 1CT 5 . 
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